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Hardware 

•NASA Advanced Supercomputing Facility enables fine-grid 
unsteady turbulent flow simulations. 
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Software 



. US3D code - Unstructured, parallel, finite-volume Navier-Stokes 
solver for thermo-chemical non-equilibrium hypersonic flows 

• Detached Eddy Simulations (or WMLES) - Hybrid RANS/ LES 
approach 


• Optional low-dissipation numerics 



Wind tunnel model (1 mm high) compared 
with the Shuttle trip (6.35 mm high) 
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X 

Unstructured hexahedral grid 
near the trip 
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Inviscid Flux Evaluation 



— + V-(F- F)= W 
dt 


Navier-Stokes equations for a 
mixture of species in 
conservation law form 


Steger-Warming Flux Vector Splitting 


F = F* +F~ = R'fCRU + R-'A-RU 


F i + 1/2 = ( R ~ l ^ R )iUi+(R~ l A~R) i+ l U M At a cell face 
Modified Steger-Warming Flux Vector Splitting 

F MI2 = {R- l KR) MI2 UMR- l A-R) MI2 U i+l Lower 

dissipation! 

A = (A±|A|)/2 
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Low Dissipation Flux 



Existing low-dissipation schemes to predict small turbulent eddies: 

- Skew-symmetric form to reduce the amplitude of the aliasing errors 

- Entropy function for secondary conservation law 
Present approach: Kinetic energy 

- Kinetic energy transport eq. derived from density and momentum eqs. 

P ? 

- Combined with internal energy: pE = pC v T + 

Kinetic energy consistent flux 

2 ( Pku'S) f = + V,.V j+1 + w t w i+i )(pu'S) f 

f f z 

Dissipative term from the modified Steger-Warming flux 




cx = 


( v 5)' 

\ ,vi • 


ICCFD6_yoon 


Implicit Time Integration 



Baseline Scheme 


Point Relaxation 

( I + Kv 2 J - A MI 2,j + Km I 2 ~ Kj-V 2 )»,J = H tJ 

+K-V 2^,-U ~ A M! 2 ,pU i+ XJ + Kj-V $U i4 - 1 - KM! ^UM 

For laminar sublayers only 

Line Relaxation 

(Modified Steger-Warming flux) 

Km: £Km + v + Ku 2 J - A lv 2J + Km: 2 - Km: 2 Wj - Km: $ U :m = H >J 




Data-Parallel Relaxation 


• Original Gauss-Seidel not parallelizable 

• Replace the G-S sweeps with a series of 
point Jacobi-like iterations or a series of line 


( A" - C " ) 8U k = H n - C n dU kA relaxations. 

QTTn+ i o TJ kmax • Second order dual time stepping (KEC flux) 

ou = ou 
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Detached Eddy Simulation 

Spalart-Allmaras 1-eq. model with compressibility correction 

°^ V = V •(- pVv)+ V • (- A /pvV A /pv)+ C " 2 V- pv ■ V A /pv + c bi Spv - c J w p( V y 
Dt o o o a A ^ 

diffusive transport production 

destruction 

• RANS approach over-predicts the turbulent dissipation levels in separated 
flow regions. 

• DES introduces a new length scale d = min(d,C DES A) 

- Model behaves like a subgrid scale model for LES away from the wall. 
-Wall-Modeled LES 
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Grid Refinement Study 





Coarse grid (15 M cells) 
solution using the modified 
S-W flux 



Medium grid (40 M cells) 
solution using the modified 
S-W flux 


Fine grid (44 M cells) 
solution using the modified 
S-W flux (temperature 
contours) 


Grid requirements for DNS: 

IB cells to resolve Kolmogorov dissipative scales 


Grid points redistributed 
(grid density 2x) 


ICCFD6_yoon 


11 



6 



Fine grid solution using KEC flux with Fine grid solution using KEC flux without a 
DES turbulence model turbulence model (ILES?) 

> The accuracy of the turbulent flow solutions is driven primarily 
by the grid resolution. 
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Formation of Upstream Vortices 




4 mm from the 
center of the trip 

(near the leading 
edge) 


3 mm from the 
center of the trip 


At the center of 
the trip 
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Conclusions 


• DES / WMLES/ Hybrid RANS-LES can be a useful tool for 
predictions of hypersonic boundary layer transition to 
turbulence triggered by an isolated roughness element. 

• It is necessary to use the low-dissipation kinetic energy 
consistent scheme on a sufficiently fine grid for an accurate 
simulation of transition. 

• Accuracy of the turbulent flow solutions is driven primarily by 
the grid resolution. 

• Interaction of vortices leads to vortex breakdown and hence 
turbulent flow. 

• Computational results agree well qualitatively with the 
experimental observations. 
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